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Transport by helical edge states of a quantum spin Hall insulator is experimentally characterized 
by a weakly temperature-dependent mean free path of a few microns and by reproducible conduc¬ 
tance oscillations, challenging proposed theoretical explanations. We consider a model where edge 
electrons experience spatially random Rashba spin-orbit coupling and couple to a magnetic impu¬ 
rity with spin S > 1/2. In a finite bias steady state, we hnd for 5 > 1/2 an impurity induced 
resistance with a temperature dependence in agreement with experiments. Since backscattering is 
elastic, interference between different scatterers possibly explains conductance fluctuations. 
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Introduction: During the last decade, the quantum 
spin Hall effect (QSHE) [1-3] has become an impor¬ 
tant example of a topologically ordered state with time- 
reversal invariance. One of its key features is the ex¬ 
istence of helical edge states [4] with a quantized con¬ 
ductance of e 2 /h per edge, as edge electrons are pro¬ 
tected from elastic single particle backscattering by time- 
reversal symmetry [1, 4, 5]. Soon after the theoretical 
prediction [6], the QSHE was realized in HgTe/CdTe 
quantum wells [7], and the quantized conductance [7] as 
well as the demonstration of nonlocal transport [8] were 
crucial signatures for this first-time experimental obser¬ 
vation. However, already in this first as well as in subse¬ 
quent experiments [7-12], deviations from the quantized 
conductance with a weak temperature dependence were 
found for edges longer than approximately 1 fj, m. More¬ 
over, in short samples, where the conductance is essen¬ 
tially quantized, small conductance fluctuations are ob¬ 
served as the back-gate voltage is tuned [7, 8, 12, 13]. Af¬ 
ter the prediction [14] of the QSHE in InAs/GaSb/AlSb 
quantum wells, the same qualitative behavior of the con¬ 
ductance as in HgTe/CdTe was observed also in these 
devices [15-18]. 

A multitude of other mechanisms beyond elastic sin¬ 
gle particle backscattering have been proposed as possi¬ 
ble explanations for the relatively short mean free path 
[19]: inelastic single particle [20-23] and two-particle 
backscattering [1, 4, 5, 23-25] which can be caused by 
electron-electron or electron-phonon interactions, both of 
which are usually considered in combination with other 
time-reversal invariant perturbations; tunneling of elec¬ 
trons into charge puddles caused by inhomogeneous dop¬ 
ing, giving rise to inelastic single particle backscatter¬ 
ing [26, 27]; coupling of edge electrons to a spin bath 
which gets dynamically polarized [28], thus effectively 
breaking time-reversal symmetry and giving rise to elas¬ 
tic backscattering in conjunction with Rashba disorder 
[29] ; time-reversal symmetry breaking by an exciton con¬ 
densate [30]; and coupling of edge electrons to a single 
Kondo impurity [4, 31, 32], to a lattice of Kondo impu¬ 


rities [33], to a single Kondo impurity in the presence of 
homogeneous Rashba spin-orbit coupling [34, 35], or to 
several Kondo impurities with random anisotropies [36]. 
Although these mechanisms are very diverse, many of 
them have in common a pronounced temperature depen¬ 
dence, usually some power law T a with positive exponent 
a for the resistance. However, only a weak temperature 
dependence has been observed experimentally, with the 
exception of a recent study using very low excitation cur¬ 
rents [37]. In fact, in some experiments, a slight increase 
of the resistance is observed when the temperature is de¬ 
creased [7, 11, 12]. With respect to the conductance fluc¬ 
tuations, only charge puddles [26, 27] as well as coherent 
scattering between several magnetic impurities with large 
spin S > 1 and uniaxial single-ion anisotropy [38] were 
considered as possible explanations. Theories that build 
on scattering from local disorder are also supported by 
recent scanning gate microscopy experiments [39], which 
identified individual scattering centers. 

In this work, we consider scattering of helical edge 
electrons from a magnetic impurity with spin S' > 1/2 
in combination with a spatially fluctuating Rashba spin- 
orbit coupling. The latter originates from a fluctuating 
electric field in the out-of-plane direction due to disorder 
in the doping layers [40-44], From a T-matrix calculation 
accounting for the combined scattering events off these 
perturbations, we derive an effective additional coupling 
to the impurity. This coupling provides a backscattering 
mechanism which is enhanced by an increased polariza¬ 
tion of the impurity with spin S > 1/2. The polarization 
of the impurity spin is determined from the steady state 
solution to a semiclassical scattering rate equation. We 
consider the linear and the nonlinear regime. Upon en¬ 
tering the nonlinear regime with the source drain voltage 
larger than temperature, the impurity gets polarized and 
the Rashba disorder induced effective coupling leads to 
an increased resistance, thus providing a possible expla¬ 
nation for the experimental results. We assume that the 
relevant Kondo temperature is exponentially suppressed 
and well below the temperature regime studied in our 
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analysis. Since the dominant contribution to backscatter- 
ing is elastic in our model, quantum interference between 
different scatterers is possible, and can give rise to con¬ 
ductance fluctuations as observed in Refs. [7, 8, 12, 13]. 

Model: The edge states are described by 

H 0 = J dx^ a (x)cj z a p(-ivM x )'S>p(x), ( 1 ) 

where tE'-t-(x) annihilates a right-moving electron, v is the 
edge velocity, and the spin quantization axis is in the 
z direction. A disordered Rashba spin-orbit coupling is 
described by [40-43] 

Hr = J dx^ a {x)al /3 {a(x),id x }^ / 3{x) , ( 2 ) 

with the correlator (a(x)a(x')) dis = VqF(x — x') and 
F( 0) = 1 [44, 45]. Being time-reversal invariant, Hr 
does not cause elastic backscattering [1], The essential 
ingredient is the coupling of electrons to a local magnetic 
impurity with spin S via 

H S = J z S z s z + J ± (S+s~ + S~s+) + J aniso (S+ + S~)s z . 

(3) 

As usual, = 5* ± iS v , and s* = (0)0-^ 4/^(0) are 

the local spin density operators of the edge electrons. 
For Janiso = 0, Hs describes a Heisenberg XXZ coupling, 
which has an axial rotation symmetry and was the sub¬ 
ject of earlier studies [4, 31-35]. In systems with axial 
symmetric Hs, the 2 component of the total spin is con¬ 
served and the dc conductance is exactly quantized [32], 
A finite Janiso breaks the axial rotation symmetry, thus 
enabling persistent backscattering in the dc limit [46]. 

In general, a coupling to the impurity spin could also 
involve terms such as S + s + +S~s~ and S z (s + + s _ ) that 
break axial rotation symmetry. However, a microscopic 
analysis, based on an isotropic sp-d exchange interaction 
of bulk electrons with the impurity, results only in the 
terms in Eq. (3), at least for HgTe/CdTe quantum wells 
[45]. Nevertheless, as we will explain, the combined pro¬ 
cesses of impurity and Rashba disorder scattering effec¬ 
tively give rise to such additional terms. 

Effective couplings from second-order processes: We 
now consider the combined scattering from the impurity 
and the Rashba disorder by using the T-matrix T(e) = 
(Hs + H r ) + (Hs + Hr)G(c)T, with G(e) = (e- Hq)- 1 . 
The T-matrix element associated with the first term of 
Eqs. (3) and (2) is 

(m; -k, -oj T z>R (ej) |to; k, a) = 
aiJ z (to | S z | to) y' f (fc ~ fe')a_fc-fc' (fc' + k)a k '-k \ 
L\Z~L ~f V hv(k — k') Hv(k' + k) J 

( 4 ) 

(cf. Fig. 1). Here, |to; k, a) denotes a product state of the 
local moment in the S z eigenstate |m) and an electron 
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FIG. 1. a) Diagrammatic representation of the second- 
order scattering process that gives rise to elastic backscat¬ 
tering without an impurity spin flip [cf. Eq. (4)]. b) The 
corresponding first-order process from the effective coupling 
H' s [cf. Eq. (5)]. The dashed line denotes the (effective) cou¬ 
pling of the impurity and electron. The dotted line with a 
cross represents the Rashba potential. 

with helical spin a and momentum hk, L the distance be¬ 
tween the left and right reservoir, ei = avhk the energy of 
the initial state |to; k , a), and a q the Fourier components 
of a(x). The energy difference between the initial and 
intermediate state in the denominator compensates the 
matrix element from Rashba disorder in the numerator. 

With regard to scattering rates, to be discussed below, 
the second-order process described in Eq. (4) can effec¬ 
tively be described as a first-order process resulting from 
the additional anisotropic coupling 

H i s = J z ,rS z (s+ + s-). (5) 

Here, J zR = 2 s Jr\J z is an effective coupling constant 
with 77 = Vo/h 2 v 2 . Analogously to the combined pro¬ 
cess described by T ZjR , there are also second-order pro¬ 
cesses where electrons scatter from the Rashba disorder 
and from the impurity by either the J± or the J an iso term 
of Hs in Eq. (3). Again, these effects can be captured by 
considering first-order processes from the effective cou¬ 
plings J±,r(S + +S~)s z and J an i so ,R(5' + s + -f S' + s _ +H.c.), 
with Jj_,r = rjj± and J an iso,R = 2 ^/fj J an i S o, respec¬ 
tively. J±.r only renormalizes the J an iso coupling already 
present in Hs- In contrast, Janiso, r, besides renormaliz¬ 
ing J l, opens an additional scattering channel via the 
S + s + + S~s~ terms. 

Scattering rates: Since we are interested in the effect 
of the impurity on dc transport, we pursue two aims: (i) 
We want to achieve a description of the impurity spin 
in a driven steady state, and (ii) we want to compute 
the transport scattering rates determining the impurity 
induced dc resistance. To achieve (i), we calculate the 
integrated rates which characterize the scattering 

of an electron from an initial state with helical spin a into 
a final state with helical spin a' , as the impurity spin is 
flipped from |m) to Im'). are obtained by weight¬ 

ing the rates for scattering from \in;k,a) to |m , ;fc , ,o' , ) 
with the probability for the initial and final state to be 
occupied or unoccupied, respectively, and then summing 
over initial and final momenta. To calculate the indi¬ 
vidual rates, we employ Fermi’s golden rule for Hg S , in¬ 
cluding all original and effective impurity couplings. For 
temperatures and voltages much smaller than the bulk 
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TABLE I. Results for the integrated scattering rates r ))))„,. 


Process 


T x (27 tAV/3/L 2 ) 

\m; a) —»• 

1 m; -a) 

jA\{m\S z \m)\ 2 r 

Kt) 

\ m +1; 4-) 

(J± + J7 n is 0i ii)| {m -1- 1| S + m) \ 2 P 

|m;t> 

1 m- 1; 4) 

TanisOjR 1 A ~ 1| S~ \m) | 2 / t 

|m; cr) —> 

\m ± 1; cr) 

(JL iso + Jl R )|(m±l|S ± |m )| 2 


excitation gap of the topological insulator, the weak mo¬ 
mentum dependence of (to'; k', a J \Hf |m; k, u) can be 
neglected. For forward scattering with <j' = cr, we then 
find 


cancel. We provide an explicit expression in the supple¬ 
mental material [45]. For the model considered here, we 
have 0 < ( < 1 and the two limiting cases have simple 
solutions: ( = 0 implies P m = S m ^s and corresponds to 
a maximally polarized local moment, while for £ = 1 the 
impurity is completely unpolarized, i.e., P m = l/(25'+l). 
Notice that (3eV = 0 implies £ = 1, because, without an 
applied transport voltage, there is no asymmetry between 
the rates for forward and backscattering [cf. Eqs. ( 6 ) and 
(7)]. The general dependence of P m on m interpolates 
between the two limiting cases: 


(1-C)C S [A 

i -c l+2S Vc; 


( 10 ) 


O m ( /3) = (m';a\ Hf \m;a) | 2 , (6) 

where the temperature dependence is due to the inte¬ 
grated occupation factors f dE f a (l — f a ) = 1//3. Here, 
ft (/f) denotes the Fermi distribution function describing 
the occupation of right-moving (left-moving) electrons 
from the left (right) reservoir. In the case of backscat¬ 
tering a' = —a, the rates involve occupation factors 
f a { 1 — f-a), leading to a voltage dependence 


r-^(AeF) 

where 


L 2 

2nh 3 v 2 [3 


| (to'; -crI Hf Im; cr) \ 2 I a (/3eV), 

(7) 


r pGpev 

r = Pj dEfM - f-a) = aPeV e ^ eV i . (8) 


At low bias voltage, (3eV < 1, ~ 1 so that the for¬ 

ward and backscattering rates have the same tempera¬ 
ture dependence. On the other hand, when fieV 1, 
the backscattering of right movers is linearly enhanced, 
while the backscattering of left movers is exponentially 
suppressed. For detailed results, see Table I. 

Master equation: We describe the state of the im¬ 
purity by a density matrix p and assume that dephas¬ 
ing from the coupling to the electron bath is sufficiently 
strong, such that we can neglect coherences and consider 
p ~ -^m \m) (m | to be diagonal in the basis of eigen¬ 

states of S z , S z \m) = m|m). We can then proceed to 
determine the steady state of the impurity spin at finite 
temperature and under an applied transport voltage V 
from the master equation 


dtPm — / P'tn r fm'mTm), 


(9) 


where T m / m = Aa'a r m'm- In a steady, state d t P m = 0, 
and we find the recursively defined solution P m - 1 = C A. ■, 
with £ = T rn _i m /T rnm _i. £ depends on /3eV and the 
(effective) coupling constants, but not on m because the 
m-dependent matrix elements of the ladder operators 


Impurity induced resistance: The backscattering prob¬ 
ability R can be related to a scattering rate 1/r by 
multiplying with the time of flight L/v. This compen¬ 
sates a dependence 1/r oc 1/L in all scattering rates, 
due to the normalization of the plane wave with a factor 
1 fy/h. In this way, both R and the edge conductance 
G = (1 — R)e 2 /h are independent of the system size 
L , as expected for a single scattering site. For R <C 1, 
R equals the impurity induced resistance normalized by 
/i/e 2 . From the effective impurity coupling we have two 
important backscattering mechanisms [48] , and hence ob¬ 
tain 


R = 


L 

v 


1 

T z , R. 



( 11 ) 


where a Fermi’s golden rule calculation yields 


1 

TL 


T 2 A- T 

j i i a. 


aniso,R 


H 2 vL 


1 

t z ,r 


t2 S 
J z, R 

H 2 vL 

m=—S 


5-1 

P±^2 I {fn + 1| S + M \ 2p m (12a) 

m=—S 

| <m| S z | m) 1 2 P m . (12b) 


Here, p± = 1—r^ m+1 /F mjm+ i accounts for the fact that 
the dc resistance is affected only by those backscattering 
events of right-moving electrons which are not compen¬ 
sated by a subsequent backscattering of a left-moving 
electron [45]. For example, in the case of vanishing 
Rashba disorder and H$ with axial rotation symmetry, 
we find Af m+1 = r m)Tn+ i, hence p± vanishes. 

Results: Although the framework that we set up so 
far does not rely on any specific assumptions about the 
couplings in H$, it is helpful to focus on the parameter 
regime J 2 J\ Jf niso for three reasons: (i) From 
a microscopic analysis we found that this regime is ex¬ 
perimentally relevant for HgTe/CdTe quantum wells [45]. 
(ii) The importance of the Rashba disorder induced effec¬ 
tive couplings with regard to the dc resistance becomes 
particularly clear in this parameter regime, (iii) A clear 
hierarchy of couplings allows one to disentangle the dis¬ 
cussion of scattering processes. The following detailed 
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FIG. 2. The normalized impurity induced backscattering 
probability R/Ro and £, both in dependence of the ratio of 
transport voltage and temperature /3eV. Ro is the backscat¬ 
tering probability for /3eV <C 1. In both cases, S = 5/2, and 
Jz/J±_ = J±/Jlniso = 10. Ro is enhanced by two orders of 
magnitude as the dimensionless strength of Rashba disorder 
g is increased from 0 to 1. 

discussion about the relevancy of Rashba disorder for the 
impurity induced resistance leads to two important re¬ 
sults: First, while R(/3eV) is a monotonically decreasing 
function without Rashba disorder, this monotonicity is 
reversed in the presence of Rashba disorder (see Fig. 2). 
Second, Ro, the backscattering probability in the limit 
(3eV -C 1, is significantly increased by Rashba disorder. 

Let us consider first the case where Rashba disorder is 
absent, i.e., rj = 0. The only nonvanishing backscat¬ 
tering rate is 1 /t_l from Eq. (12a), which is small as 
it arises from an interplay of scattering due to J± and 
J aniso■ In particular, Rq is found to be proportional to 
the harmonic mean of J]_ and 2Jf niso , because p± = 
[1 + JI/2JLJ- 1 for j3eV = 0. Consequently, when 
J\ and Jf niso are very different in magnitude, it is the 
smaller of the two which determines the magnitude of 
Rq. With increasing /3eV, the backscattering of right 
movers via J±S + s~ becomes increasingly dominant rel¬ 
ative to other scattering rates and tends to polarize the 
impurity, such that ( ~ 2J 2 niso /fdeVJ^ approaches zero 
in the large {3eV limit (cf. Fig. 2). With increasing po¬ 
larization, the probability for the impurity spin to be in 
state IS 1 ) increases, and the probability for an individual 
right-moving electron to be backscattered is suppressed. 
This leads to the monotonic decrease of R(/3eV) shown 
in Fig. 2, with R ~ </ 2 niso //3eF in the large j3eV limit 
[49]. 

Rashba disorder, described by a finite rj, has a pro¬ 
found effect on the impurity induced resistance. Since 
J 2 » J]_ » J 2 niso , l/r-,R dominates 1 /tj_ already for 
very weak Rashba disorder with g > 4Jf niso /jf. The 
magnitude of Ro is then determined by J| R instead of 
J 2 ni SO ) which, depending on the precise values of the 
couplings and g, can be a large difference. Regarding 
the dependence on the polarization, 1/t Z) r is qualita- 


FIG. 3. A right-moving wave (gray) is phase coherently 
backscattered from impurities at Xo (red) and Xo + l (black). 
Destructive (constructive) interference occurs when IcfI = vn 
with half-integer (integer) v\ cf. a) [b)]. 


tively different from 1/tj_. Evaluating the sum over m 
in Eq. (12b) for the limiting cases of perfect polarization, 
Pm = S mj s, and vanishing polarization, P m = l/(25 + l), 
yields S 2 in the former and S(S + l)/3 in the latter 
case, respectively. This shows that, for spin S > 1/2, 
the rate 1 /t Zj r increases with increasing polarization, be¬ 
cause S 2 > S(S + l)/3. However, l/r z r is independent 
of £ for S = 1/2. Thus, in contrast to the case with¬ 
out Rashba disorder, R is now found to monotonically 
increase with (3eV when S > 1/2. 

For real samples with several impurities, the total 
backscattering rate is proportional to the number of im¬ 
purities in the absence of localization. Based on Eq. (11), 
we estimate the mean free path to be 4 pm in 7.0 nm wide 
HgTe quantum wells with a lattice constant a = 0.65 nm 
and v = 4 x 10 5 m/s, by assuming (i) that there is a con¬ 
centration of 10~ 4 Mn 2+ ions (S = 5/2) per unit cell, (ii) 
£ = 0.5 and spatial average (J 2 /h 2 v 2 ) = 0.035 [45], and 
(iii) rj « 3 based on the estimate for Vo from Refs. [29, 44]. 

A plane wave e zkFX with Fermi momentum fcp can 
be backscattered elastically from impurities at Xq and 
xo + l by the process Eq. (4) (see Fig. 3). The partial 
waves e - ifc F(x-x°)+ifepx o anc j e -ik F (x-x 0 -l)+ik ¥ ’(x 0 +l) j n _ 

terfere, resulting in a contribution to the backscattering 
probability oc cos kyl. Tuning fcp via a back-gate volt¬ 
age alternately causes constructive or destructive inter¬ 
ference. Hence, conductance fluctuations occur, in agree¬ 
ment with experiment. Several impurities give rise to 
more complex oscillations. 

Summary: We determined the dc resistance of heli¬ 
cal edge states in the presence of Rashba disorder and a 
magnetic impurity with spin S > 1/2. As a key result, 
we find that combined scattering from both Rashba dis¬ 
order and impurity is enhanced as the impurity becomes 
more polarized, giving rise to a resistance that slowly 
increases as the ratio of transport voltage and tempera¬ 
ture increases, in agreement with experiments [7, 11, 12]. 
Since backscattering is elastic, quantum interference can 
explain the occurrence of conductance fluctuations. 
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SUPPLEMENTAL MATERIAL 


HgTe/CdTe quantum wells: microscopic considerations about Hs 


In many cases, it is a useful simplification to consider the edge states of a topological insulator as one-dimensional 
objects by neglecting the spatial extension of their wave function into the bulk. Here, we take this spatial extension 
as well as the orbital character of the edge states into account, in order to derive from the exchange interaction of 
bulk electrons in HgTe or CdTe with a point-like magnetic impurity, the effective exchange interaction between such 
an impurity and electrons in the helical edge states of a topologically non-trivial HgTe/CdTe quantum well. In what 
follows, we choose coordinates such that the x-axis coincides with the edge and the z-axis is parallel to the growth 
direction of the quantum well. The location of the impurity is chosen to be r Q = (0, yo, z 0 ). 

In Ref. [1], the BHZ Hamiltonian was derived, to describe the low energy physics of HgTe/CdTe quantum 
wells. This effective Hamiltonian for the quantum well acts on the quantum well subbands |£11,±) and |iJl,±), 
which are the subbands with the lowest energy at the T-point. Hence, it is usually written as a 4 x 4 matrix 
in the basis (Iffl, +), \H1, +), \E1, —), \H1, —)). The BHZ Hamiltonian was derived from the six-band bulk Kane 
model, which is formulated in the basis (|T 6 ,+1/2), |T 6 , —1/2), |T 8 ,+3/2), |r 8 , +1/2), |T 8 , —1/2), |T 8 , —3/2)). As 
a consequence |£T,±) and |ffl,±) can each be represented in the basis of the six-band Kane model: |i?l,+) = 
/i(z) |T 6 , +1/2) + / 4 (z) |T 8 , +1/2), \H1, +) = t /* 3 (z) |T 8 , +3/2) [ 1 ]. Edge states, which can be found as eigenstates of 
the BHZ Hamiltonian with appropriate boundary conditions [2, 3], thus can be expressed in the basis of the six-band 
Kane model. To be precise, we express the edge state | k,a) in position space 

(r | k,a) = ~^=e lkx Mbuik{z)M qv ,(k, y) \a) , (13) 

where the spinors |f) = (1,0) T and |j.) = (0,1) T correspond to right- and left-movers, respectively. The matrices 
Mbuik and M qw explicitly read 



//lO) 

0 

0 

0 ^ 





0 

0 

-/?(*) 

0 


^ cei 

°\ 

Afbulk (z) = 

0 

h{z) 

h{z) 

0 

0 

0 

0 

0 

, M qw (k,y) = g k (y) 

chi 

0 

0 

C E1 


0 

0 

/!(*) 

0 


V 0 

chi) 


V 0 0 0 -/ 3 *(z)/ 


The 4x2 matrix M qw is used to express the spinors |er) in the basis of the BHZ Hamiltonian, while the 6x4 matrix 
Afbuik expresses the elements of this basis in the basis of the six-band Kane model. The components cei and chi can 
be obtained analytically [2, 4], we use c E1 = 0.13 and c'j 11 = 0.87 independent of k. In addition, gk{y) describes the 
exponential decay of the edge state wave function away from the edge; it is known analytically [2, 3]. The functions 
fi(z ) were described in Ref. [1]. fi and /3 are even and real, whereas /4 is odd and imaginary; all three functions 
decay exponentially in the CdTe region. 

The interaction of bulk electrons in CdTe and HgTe with the local moment of a magnetic impurity like a Mn 2+ ion 
may be described by the Hamiltonian = aS • sr 6 + (/3/3)S • Ir 8 , [5]. Here, S acts on the local moment, sr 6 is 

the local spin density of electrons in the Tg band at the impurity site, Jr 8 is the density of total angular momentum 
of holes in the T s band evaluated at the impurity site, and a and /? are exchange constants. Expressing in the 

basis of the six-band Kane model, |T 6 , +1/2), |T 6 , —1/2), |T 8 , +3/2) , |T 8 , +1/2), |T 8 , —1/2), |T 8 , —3/2), and keeping 
in mind that the operators S + , S~ , and S z act on the impurity state, we find 


H. 


bulk 



V 



£1 

3 2 


/ 3 s z 

a /3 S+ 
0 

V 0 


V3S~ 

s z 

2 S+ 

0 


\ 


0 

2 S~ 

- s z 
735 ’+ 


0 \ 

0 

V3S~ 
~3S Z J ) 


6{x)S{y - y 0 )S(z - z 0 ). 


(15) 


We are interested in the effective exchange interaction H^ se between the edge electrons and the impurity. To obtain 

bulk M bu ik to 


the corresponding 2x2 Hamiltonian in the basis (|t), |j_)), we first calculate H^(x,y) = f dz M^ ulk ET/ x 
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find the effective exchange interaction in the basis of the BHZ Hamiltonian and afterwards compute H^ ge (x) = 
JdyM;[ w H^M qw . 


with 


H^(x,y) = J dzMl ulk H^ lk M hulk 


(JTS‘ (JD*s+ 

j-qw g— jqw g Z 

J 3 qw S+ 0 

V 0 0 


{j$ w ys~ 

o 

-jrs* 

-(J 4 qw )*S+ 


0 \ 
0 


-J 2 w s~ 

-J2 w s z / 


S(x)S(y-y 0 ), 


J i W ( z o) = ^(3a|/i(2o)| 2 + P\h(z 0 )\ 2 ), 
JT(zo) = ^|/ 3 (-o)| 2 , 

J 3 W (~o) = ^(-3a/i(2 0 ) 2 + 2/3/ 4 (2 0 ) 2 ), 

JTM = ^PfZ( z o)U(z 0 ). 


^e e x dge W= / dy M qw 


JzS~ + Janiso'S' + Ja n iso'-’ + 

J| 5+ 


J* ± S~ 


~(J Z S~ + Janiso S 1 + Janiso) , 


<50*0 


with 


Jz(VO,ZO ) = |5fc(yo)| 2 (^i W |cei | 2 + J£ W |c/n| 2 ) 

Jl_ (2/0,2o) = l5fc(j/o)| 2 -/3 Wc ll 

Janiso (2/0 j A)) = l<7fc(j/o)| J 4 hElCffi 


(16) 


(17a) 

(17b) 

(17c) 

(17d) 

(18) 

(19a) 

(19b) 

(19c) 


is equivalent to defined in Eq. (3) of the main text, except that the complex phase of J an i so was neglected 
in the main text for brevity; J± is real. Values for a and (3 can be found in the literature for HgTe as well as CdTe 
[5, 6]. J z is found to be the largest of the three couplings independent of zq. fs(z) is odd such that J an iso vanishes if 
the impurity is located at zq = 0, i.e. in the middle of the quantum well. Consequently J? S> Jj_ Jf niso is realistic 
for Zo ~ 0 . 

The reason why a finite Janiso can arise from an isotropic exchange interaction like lk can be understood from 
the calculation and was explicitly stated in Ref. [7]: The electrons in the edge states in HgTe/CdTe quantum wells 
are not spin-polarized because of the spin-orbit interaction in these materials. This also holds for InAs/GaSb/AlSb 
quantum wells. 


Details and intermediate steps regarding the effective coupling J z ,r 


In this section we provide the details of the determination of the effective coupling J z ,r in Eq. (5) of the main text, by 
calculating the disorder averaged squared modulus of matrix elements of the T-matrix term T z ^(e) = H R G(e) J z S z s z + 
J z S z s z G(c)Hr for elastic backscattering. Let us calculate explicitly the matrix element of the first term, which is 
represented in Fig. 1 a) of the main text, the contribution from the second term is obtained analogously. We use the 
Fourier representations J z S z s z = £ S z T, kl ,k 2 and H R = E fel ,fc 2 ( fc i + fe 2H l -fc 2 4 1 ,a f7 a l S Cfc =’' 8 to find 

(m; —k, —a\ Ur- 


CTvhk — Hq 


s \m:K,cr) 


= E 

k' 

= £ 


(m;—fc ,—+ k 2 ) a kl -k 2 c { ua Ck 2 ,/}\m-, k',a)(m; k', <r\^S z a z aoi T, kl ,k 2 c\ ua c k ^ a \m-,k,v) 

crhv(k — k ') 

(m; -k, -a\ ^ <r (-fc + fc , )a_ fe _ fc /cL fe) _ 0 .c fc q 0 .|m; k',a)(m;k',a\dg-S z a z a cl, a c kk ,cr\m-, k,cr) 


ahv{k — k') 




( 20 ) 












The second term yields almost the same result, except that the sum now is ak'-k ■ Hence for the disorder averaged 
squared modulus of the matrix element, which enters into the scattering rates, we find 


<| (m; -k, cr| T 2 , r |m;fc,cr) | 1 2 ) dis = (H S z \m) | 2 (^(a fc /_ fc + a_ fc _ fc ,) y^(a fc /_ fc + a_ fc _ fc /)*) dis 


L 3 4 5 6 7 fi, 2 1! 2 

= Y2^ r lJz \ ( m l S Z \ m ) | 2 


( 21 ) 


where rj = Vo/h 2 v 2 . We used (a g a*,) d i s = Va5 q ^ q 'F(q), with F(q ) = j- Ax F(x)e lqx and F( 0) = 1. Since 
| (m; — k, —cr| J 2j rS' 2 (s + + s~) \m;k,a) | 2 = J 2 R | (to| S' 21 |to) | 2 , we can capture the effects from the second order 
process on the scattering rates by considering first order processes from the effective coupling Eq. (5). 


Explicit expressions for £ and p±. 

As mentioned in the main text £ = T m-i,m Since | (to — 1| S~ |to) | 2 = | (to| S + | m — 1) | 2 , these matrix 

elements cancel. Consequently £ does not depend on m, but only on the (effective) coupling constants and on /3eV 
through the integrals I a 


2 (^niso + J±,r) + + ■^mso.r) + -^anisn.R 

2 (^aniso + ^1 ,r) + ^^aniso.R + ^ (J± + ^fniso.R) 


The probability p± is calculated according to the following considerations: Assume a right-moving electron is 
backscattered while flipping the impurity from m to to + 1, i.e. assume that a backscattering event of the type 
|m;t) —> |to + 1;£) occurred. Then, because the impurity is assumed to be in a stationary state, at some time 

another scattering event |to + 1; cr') —» |to; cr) has to take place. If the second event happens to be a backscattering 
of a left-moving electron (cr' =4_, and cr =f), then in total, the two scattering events left the current transmitted 
between the two reservoirs unchanged. However, if the second event is a forward scattering event (cr' = cr) or 
the backscattering of another right-mover (cr' =f, and cr =|), then the first backscattering event indeed had a net 
effect on the current and thus contributed to the dc resistance. Since it does not matter which momentum the 
states involved in the scattering have, the probabilities for the two different cases are proportional to the integrated 
scattering rates T™ m ,. The probability that the first scattering event does not have a net effect on the resistance is 
Pm+i m a m+v < - )n the °ther hand, the probability that the first scattering event does have a net effect on the 
resistance is ?C+ 1|TO a T^ m+1 + T^ m+1 + T^ m+1 . Since fC+ ljm + pZh, m = 1 we find 


DOS 

Pm+1- 


1 + 


r 


itt 

m,m+1 


pti 

m,m+1 

I "p-i-i- p-i-T 

' 1 m,m+1 ' L m,m+1 


-1 


(23) 


Since r^ CT m+1 oc | (to| S \m + 1) | 2 independently of cr and cr', m is independent of to- and we can define 


P -l 




(■ Jl+Jl 


aniso,R 


)/•*■ 


2 (^anis 


J 1,r) + ^aniso.R^ 


-1 


(24) 


[1] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, 314, 1757 (2006). 

[2] B. Zhou, H.-Z. Lu, R.-L. Chu, S.-Q. Shen, and Q. Niu, Phys. Rev. Lett. 101, 246807 (2008). 

[3] M. Wada, S. Murakami, F. Freimuth, and G. Bihlmayer, Rev. B 83, 121310 (2011). 

[4] P. Virtanen and P. Recher, Rev. B 85. 035310 (2012). 

[5] I. A. Merkulov and A. V. Rodina, Exchange Interaction Between Carriers and Magnetic Ions in Quantum Size Heterostruc¬ 
tures, In: J. A. Gaj, and J. Kossut eds., Introduction to the Physics of Diluted Magnetic Semiconductors (Springer, Berlin, 
Heidelberg, 2010). 

[6] J. K. Furdyna, Appl. Pins. 64, R29 (1988). 

[7] J. Maciejko, T. L. Hughes, and S.-C. Zhang, Annual Review of Condensed Matter Physics 2, 31 (2011). 







